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ABSTRACT 

In an attempt to probe the magnetic field morphology near the massive young star Orion-IRc2, we 
mapped the linear polarization of its J=2-l SiO masers, in both the v=0 and v=l vibrational levels, 
with 0.5" resolution. The intense v=l masers are confined to a narrow zone 40 AU from the star. Their 
polarization position angles vary significantly on time scales of years. For the v=l masers the stimulated 
emission rate R is likely to exceed the Zeeman splitting gVL due to any plausible magnetic field; in this case 
the maser polarization need not correlate with the field direction. The much weaker v=0 masers in the 
ground vibrational level lie 100-700 AU from IRc2, in what appears to be a flared disk. Their fractional 
polarizations are as high as 50%. The polarization position angles vary little across the line profile or 
the emission region, and appear to be stable in time. The position angle, PA. = 80°, we measure for 
the J=2-l masers differs by 70° from that measured for the J=l-0 SiO transition, possibly because of 
Faraday rotation in the foreground, Orion A, H II region. A rotation measure RM = 3.3 x 10 4 rad m~ 2 
is required to bring the J=2— 1 and J=l-0 position angles into concordance. The intrinsic polarization 
position angle for both transitions is then 57° , parallel to the plane of the putative disk. Probably the 
magnetic field threads the disk poloidally. There is little evidence for a pinched or twisted field near the 
star. 

masers — polarization — stars: formation — stars: individual 



Subject headings: magnetic fields 
(Orion-IRc2) 

1. INTRODUCTION 

As an interstellar cloud collapses to form a star, the 
magnetic field lines threading the cloud are expected to 
be dragged along with the infalling gas, producing an 
hourglass-shaped field pattern centered on the star (e.g., 
Galli & Shu 1993). In many theoretical models the field is 
then wound up by rotation in a circumstellar disk, produc- 
ing a toroidal field component which collimates the out- 
flow from the star into two oppositely directed lobes; the 
outflow carries away angular momentum, allowing further 
accretion onto the star (Uchida & Shibata 1985; Pudritz 
& Norman 1986; Newman, Newman, & Lovelace 1992). 

So far the observational evidence for pinched or twisted 
fields near young stars is scant, although there are possible 
detections toward W3 (Greaves, Murray, & Holland 1994), 
Mon R2 (Greaves, Holland, & Murray 1995), and a num- 
ber of other cloud cores (Greaves & Holland 1998). These 
observations used the linearly polarized emission from cir- 
cumstellar dust to probe the magnetic field morphology. 
Spinning dust grains tend to align their long axes perpen- 
dicular to the magnetic field, so thermal emission from the 
grains is polarized perpendicular to the field. 

Toward the Orion molecular cloud, Schleuning (1998) 
found evidence for a large scale hourglass distortion of 
the magnetic field from dust polarization observations at 
100 and 350 /im. The pinched region has a radius of 
~0.5 pc (200") and is centered 0.1 pc southwest of the 
Kleinmann-Low nebula, a site of high mass star forma- 



tion. On smaller scales, within the Orion-KL complex, 
Chrysostomou, et al. (1994) found evidence for a twisted 
field within roughly 1500 AU (3") of the star IRc2, based 
on imaging polarimetry of the 2 ^jm S(l) line of H2. The 
line emission, assumed unpolarized, originates from shock- 
excited H2 in the bipolar outflow from IRc2. Absorption 
by aligned grains in front of the outflow produces the po- 
larization. Aitken et al. (1997) also found evidence for 
a twisted field near IRc2 from 1.5" resolution polarization 
maps at 12.5 ^m and 17 /im. Spectropolarimetry of the 8- 
13 and 16-22 /im bands was used to separate the effects of 
absorption and emission, which lead to polarizations par- 
allel and perpendicular to the magnetic field, respectively. 
Rao, et al. (1998) obtained 4" resolution interferometric 
maps of Orion-KL at wavelengths of 1.3 and 3.5 mm, where 
the dust emission is certainly optically thin. These maps 
showed an abrupt 90° change in the polarization angle 
southeast of IRc2, in agreement with the infrared data. 
But because the region of anomalous polarization was sus- 
piciously coincident with the redshifted lobe of the bipolar 
outflow from IRc2, Rao et al. suggested that dust grains 
in this zone might be aligned by a wind from the star, 
rather than by the usual Davis-Greenstein mechanism. If 
the grains are mechanically aligned then thermal emission 
can be polarized parallel to the magnetic field, so there is 
no need to conclude that the field is kinked. 

Unfortunately, polarized dust emission from IRc2 may 
not be a reliable probe of the magnetic field near the cen- 
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tral star because of confusion by dust emission from the 
surrounding molecular cloud. The stellar position is coin- 
cident with radio source T (Menten & Reid 1995). The 
IRc2 infrared peak is offset <~1" northwest of source I, 
while the millimeter- wavelength peak is offset south- 
east of source I, toward the 'hot core' clump. The radio 
spectrum of source I itself from 8 GHz to 86 GHz is consis- 
tent with free-free emission from ionized gas in a circum- 
stellar shell; there is no evidence for excess flux at 86 GHz 
which would be the signature of circumstellar dust (Plam- 
beck et al. 1995). 

By contrast, SiO masers do unambiguously originate in 
gas very close to source I. Intense masers in the v=l vi- 
brational level are confined to a zone approximately 40 
AU from the star (Plambeck, Wright, & Carlstrom 1990; 
Menten & Reid 1995; Greenhill et al. 1998; Doeleman, 
Lonsdale, & Pclkey 1999), while weaker masers in the 
ground vibrational state are found 100-700 AU from the 
star, in an elongated, bowtie-shaped region (Wright et al. 
1995; Chandler & de Pree 1995). Both the v=l and v=0 
masers are linearly polarized (Barvainis 1984; Tsuboi et al. 
1996). Since the maser polarization direction is expected, 
under many circumstances, to be parallel or perpendicular 
to the magnetic field direction (Goldreich, Keeley, & Kwan 
1973), SiO masers may provide the best probe of the mag- 
netic field morphology within 700 AU of source I. Indeed, 
from early measurements showing a symmetric rotation of 
the polarization angle across the double-peaked v=l maser 
spectrum, Barvainis (1984) concluded that the magnetic 
field direction rotated with azimuth in a disk around the 
star. 

Here we discuss polarization measurements of both the 
v=0 and v=l SiO masers toward Orion- IRc2 obtained with 
0.5" angular resolution, ^250 AU at the distance of Orion. 
We find that the polarization of the v=l maser is time 
variable, casting doubt upon its usefulness as a probe of 
the magnetic field. However, the v=0 maser polarization 
appears to be a good probe of the field morphology. It 
provides little evidence for pinched or twisted fields near 
the star. 

2. OBSERVATIONS 

Observations were made between 1996 June and 2001 
March with the BIMA array. 1 Orion-KL was observed 
using the A, B, and C antenna configurations, providing 
antenna separations from 6.3 m to 1.8 km. The correla- 
tor was configured to allow simultaneous observations of 
the J=2-l SiO transitions in both the v=l and v=0 vi- 
brational states, at frequencies of 86.243 and 86.847 GHz. 
The velocity resolution was 0.34 km s _1 . Single sideband 
system temperatures, scaled to outside the atmosphere, 
ranged from 170 K to 350 K. The data were calibrated 
and analysed using the MIRIAD software package (Sault, 
Teuben, & Wright 1995). 

The receivers are sensitive to a single linear polariza- 
tion, but are equipped with movable quarter wave plates 
to observe right or left circular polarization. As described 
by Rao, et al. (1998), each receiver was switched between 
LCP and RCP with a Walsh function pattern in order to 
sample all possible crosscorrelations (RR, RL, LR, LL) on 
each baseline as rapidly as possible. The integration time 

1 Operated by the Berkeley-Illinois-Maryland Association with supp 



for each observation was 11.5 seconds. The data were self- 
calibrated on the strongest v=l maser feature, using the 
LL and RR crosscorrelations, then averaged in 5 minute 
blocks to produce quasi-simultaneous dual polarization ob- 
servations. 

The quarter wave plates are designed for a frequency 
of 89 GHz, where they have a leakage to the unwanted po- 
larization of around 1%. At 86 GHz the leakage is about 
3%. We calibrate the leakage by observing a point source 
over a wide range of hour angle. For altitude-azimuth an- 
tennas the phase of a linearly polarized source rotates with 
the parallactic angle, whereas the leakage terms are con- 
stant. We used observations of 3c273, 3c279, and the SiO 
v=l maser itself to calibrate the leakage. The measured 
polarizations of both the v=l and v=0 masers are quite 
insensitive to the leakage correction: for the weakly po- 
larized v=l line, the source is unresolved and the leakage 
averages out over the wide range of parallactic angles ob- 
served, while for the v— line the source polarization is 
much larger than the leakage. 

Images were made for a range of weightings of the uv- 
data. With uniform weighting the synthesised beam is 
0.8" x 0.3" FWHM. For the final images we used robust 
weighting, deconvolved using the CLEAN algorithm, and 
convolved with a circular 0.5" FWHM Gaussian beam. 

3. RESULTS 

3.1. The v=l Masers 

Figure 1 summarizes the observational results for the 
v=l masers over a 3.5 year period. For each of 6 epochs, 
we plot the total intensity, fractional polarization, and po- 
larization position angle across the double-peaked spec- 
trum. All 3 quantities vary on time scales of years. The 
fractional polarization usually is less than 10%, and tends 
to be anticorrelated with maser intensity - that is, the 
strongest maser features typically have the lowest frac- 
tional polarization. The weak maser feature near Vlsr 7 
km s -1 which appeared in 1999 has an anomalously high 
fractional polarization of 30%. The polarization position 
angles vary considerably across the line profiles, although 
a crude median value is P. A. ~ 80°, shown by the dashed 
line in Figure 1. Often the position angle "kinks" at the 
peak of a maser feature. The symmetric pattern of po- 
larization position angles seen in the 1996 August data 
is similar in character to that measured by Barvainis & 
Predmore (1985) in 1981 June. 

The v=l masers are unresolved in 0.5" resolution im- 
ages. Nevertheless, the high signal to noise ratio of the 
maps, in excess of 5000:1 on the strongest maser features, 
allows us to fit the centroid of the maser position in each 
velocity channel to an accuracy of ±0.01". The centroids 
are arranged in an 0.15" diameter ring with a stable, 
systematic velocity pattern (Plambeck, Wright, & Carl- 
strom 1990; Wright et al. 1995; Baudry, Herpin, & Lucas 
1998). In Figure 2, the polarization vectors for individual 
1 km s _1 velocity channels are plotted on the correspond- 
ing centroid positions for 3 different epochs. Variations 
in the polarization angle do not appear to correlate with 
movements of spots. Below we argue that the polariza- 
tion angles of the v=l masers may be determined by the 
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maser beaming direction, rather than the magnetic field, 
accounting for the time variability. 

Plambeck, Wright, & Carlstrom (1990) were able to fit 
the v=l maser spectrum and centroid positions with a 
model of maser emission from a rotating, expanding disk. 
Higher resolution, VLBA observations (Grccnhill ct al. 
1998; Doeleman, Lonsdale, & Pelkey 1999) show, however, 
that the strongest maser features are clustered into four 
groups, plausibly in the walls of the bipolar outflow from 
source I. Although Doeleman et al. emphasize the dif- 
ferences between the disk and outflow models, the distinc- 
tions are in part semantic - the disk modeled by Plambeck 
et al. had a sharp inner radius of 40 AU and a thickness of 
40 AU, and the maser pump rate and SiO density were fit- 
ted to obtain particularly strong maser emission from the 
inner walls of this cylinder, not so different from the out- 
flow model. Nevertheless, it is clear that the velocity field 
hypothesized by Plambeck et al. is an oversimplification. 

3.2. The v=0 Masers 

The bulk of the v=0 J=2-l SiO emission from the 
Orion-KL core is thermal. Much of it originates from 
high velocity gas in the bipolar outflow from IRc2 (Wright 
et al. 1983). With arcsecond angular resolution, most of 
this extended emission is resolved out, revealing a bright 
"bowtie" centered on source I. The peak brightness tem- 
peratures in the bowtie are > 2000 K, almost certainly 
indicating maser emission (Wright et al. 1995). Figure 3 
compares total intensity and polarized spectra of Orion- 
IRc2 generated from high and low resolution maps. The 
polarized flux densities integrated over 3" x 3" and 30" x 30" 
boxes are nearly identical, indicating that most of the po- 
larized emission originates from the masers. 

SiO channel maps with 0.5" resolution are shown in Fig- 
ure . The r.m.s. noise in the Q and U Stokes intensity 
maps is ~0.05 Jy/beam. Polarization vectors are over- 
laid wherever the linearly polarized flux P = (U 2 + Q 2 ) 1 ^ 2 
is detected with 4ct or greater significance; hence the po- 
larization position angles shown are uncertain by A(f> = 
AP/(2P) < 7°. Apart from a slight twist of the polariza- 
tion vectors at the northern edge of the Vlsr 12.5 km s _1 
map, there arc no obvious systematic changes in polar- 
ization angle with velocity or position. The polarization 
direction, P. A. ~ 80°, is offset by 30° from the long axis of 
the bowtie. The fractional polarization is as high as 50% 
in the channels near the upper and lower edges of the line 
profile, near Vlsr —5 and +15 km s _1 . 

The SiO linewidths in the bowtie are too large for this 
gas to be gravitationally bound to source I. The "shell" 
H2O masers in Orion-IRc2 overlie the inner regions of 
the bowtie. Proper motion measurements show that these 
masers are moving outward from source I at ^20 km s , 
perpendicular to the axis of the high velocity outflow 
(Greenhill et al. 1998). Probably the v=0 SiO masers and 
H 2 masers are part of a low velocity outflow which is 
pushing into, or around, dense gas in the equatorial plane 
of the star, plausibly the remnants of the "pseudodisk" 
which formed in the intial collapse of the cloud (Galli & 
Shu 1993). 

4. DISCUSSION 
4.1. Maser Polarization Mechanisms 



It is remarkable that a magnetic field with a strength 
of a few milligauss can affect an SiO maser in any ob- 
servable way, given that the Zeeman splitting gVL in- 
duced by the field is orders of magnitude smaller than the 
maser linewidth Av. For SiO, gVL/2-K ~ 0.2 Hz/milligauss 
(Nedoluha & Watson 1990), while individual J=2-l maser 
features have linewidths of order 0.5 km s _1 , or 1.4 x 10 5 
Hz. Nevertheless, polarization arises because of the selec- 
tion rules which govern dipole transitions between energy 
levels. Radiation polarized parallel to the field cannot 
change the component of molecular angular momentum 
m along the field direction, hence Am = 0; conversely, 
Am = ±1 for radiation polarized perpendicular to the 
field (Townes & Schawlow 1975). These rules cause the 
maser gain to differ for polarizations parallel and perpen- 
dicular to the magnetic field. 

A simple example given by Goldreich, Keeley, & Kwan 
(1973) illustrates this mechanism. Consider a maser in the 
J=l transition. In the presence of a magnetic field, the 
J=l state is split into 3 sublevels (m = +1,0,-1). If a 
maser beam propagates at right angles to the field, radia- 
tion polarized parallel to the field induces transitions only 
from the ( J, m) = (1,0) state to the (0,0) state, while ra- 
diation polarized perpendicular to the field induces transi- 
tions from both the (1, +1) and (1,-1) states to the (0, 0) 
state. If the pumping rates into the 3 upper sublevels are 
equal, and if the maser is saturated - that is, every exci- 
tation to the upper state leads to the emission of a maser 
photon - then two photons polarized perpendicular to the 
field are emitted for every one photon polarized parallel 
to the field. Hence the net fractional linear polarization is 
1/3, perpendicular to the field. 

In the more general situation where the maser prop- 
agates at angle < 9 < ir/2 with respect to the field, 
radiation polarized parallel to the projected field direction 
couples to all 3 magnetic sublevels, while that polarized 
perpendicular to the field couples only to the m = ±1 
sublevels. Hence, as 9 decreases from ir/2 toward 0, the 
gain of the parallel component increases while the gain 
of the perpendicular component decreases. In the limit 
gQ 3> R ^ r, where R is the stimulated emission rate and 
r is the decay rate from collisions or spontaneous emission, 
Goldreich, Keeley, & Kwan (1973) show that the net lin- 
ear polarization for a J=l-0 maser is perpendicular to the 
projected field for 9 > 54.7° (sin 2 9 > 2/3), and parallel to 
the projected field for 9 < 54.7°. 

Western & Watson (1984) found that these same asymp- 
totic limits applied to the J=2-l rotational transition. In 
addition, they considered the case where the magnetic sub- 
levels were pumped unequally, which can easily occur if the 
vibrational energy levels are excited by an anisotropic radi- 
ation field, for example by infrared photons from a nearby 
star. Small differences in the absolute populations of these 
sublevels correspond to large fractional differences in the 
population inversion, and hence the maser gain. There- 
fore, anisotropic pumping tends to enhance the fractional 
polarization, and also changes the angle 9 at which the 
polarization flips from perpendicular to parallel. Nonethe- 
less, as long as gQ 3> R and T, the net linear polarization 
is either parallel or perpendicular to the projected field 
direction. 

Unhappily, the requirement gtt ^> R is likely to be vio- 
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lated for intense mascrs, for which the stimulated emission 
rate R becomes large. As R increases, the "good" quan- 
tum axis shifts from the magnetic field direction to the 
propagation direction of the maser beam. Calculations by 
Ncdoluha & Watson (1990) show that the position angle 
of the linear polarization is neither parallel nor perpendic- 
ular to the projected magnetic field in this case, over a 
broad range of maser intensities. For example, Figure 3 of 
Ncdoluha & Watson shows that the polarization position 
angle for a J=2-l SiO maser propagating at 15° to the 
magnetic field direction shifts from 20° for R = gQ,, to 80° 
for R = lOO.gfi. As the maser intensity increases and R 
exceeds gfl, the calculations also show that the fractional 
polarization decreases. 

4.2. The v=l Masers 

The v=l SiO masers near IRc2 are sufficiently intense 
that the stimulated emission rate R almost certainly is 
larger than the Zeeman splitting gQ,. The brighter J=l-0 
maser spots in the VLB A map of Doeleman, Lonsdale, & 
Pelkey (1999) have flux densities of 1 Jy in an 0.2 milliarc- 
second beam, corresponding to brightness temperatures of 
2 x 10 10 K. To estimate the stimulated emission rate R for 
an SiO molecule near the surface of such a maser, one must 
guess the solid angle dfl into which the maser radiation is 
beamed. Taking dSl = 10~ 2 steradians, 



R = BU = 



8ttV (•/ + !) 



dCt 



3 * ( 2j +3 r c - 30s " (1) 

where B is the Einstein B coefficient for the J+l — > J tran- 
sition, U is the radiation energy density in ergs cm -3 Hz -1 , 
fi = 3.1 x 10~ 18 esu-cm is the dipole moment for SiO, and 
I = 2kT/X 2 is the maser's specific intensity. For compari- 
son, the spontaneous radiative decay rate from the v=l to 
the v=0 vibrational level is 5 s _1 (Hcdclund & Lambert 
1972) and the SiO-H 2 collision rate in 1000 K gas with 
density 10 9 cm~ 3 is ~1 s _1 , using the total rate coeffi- 
cient U 1.4 x 10- 9 (T/T ) - 39 cm 3 s" 1 (T = 2000 K) 
given by Bieniek & Green (1983). Since R » T, the maser 
is saturated. 

In order to fulfill the requirement gtt 3> R, magnetic 
fields >300 mG would be required. For a 300 mG field the 
magnetic pressure B 2 /8tt is 25 times greater than the ther- 
mal pressure nkT, which seems unlikely. Thus, we reach 
the unwelcome conclusion that the polarization direction 
of the v=l mascrs is a poor indicator of the magnetic field 
direction near IRc2. The model calculations of Nedoluha 
& Watson (1990) show that the polarization position angle 
is a function of R/gQ, and hence of the maser flux density. 
Thus, one can understand both (1) the variations in po- 
larization angle which occur across the maser line profile, 
and (2) the changes in polarization direction which occur 
on time scales of years as the maser intensities change. 
For R/gfl > 1 the models also suggest that the most in- 
tense masers will have the lowest fractional polarization, 
as observed. 

An observational result of note is the maser feature with 
30% fractional polarization which appeared at Vlsr~ 7 
km s _1 in 1999 June. Such high fractional polarization 
is difficult to achieve unless the masers are anisotropi- 
cally pumped by infrared radiation from the central star 
(Western & Watson 1984). Anisotropic pumping also can 
change the critical angle at which the polarization becomes 



parallel or perpendicular to the field, perhaps explain- 
ing the anomalous position angle observed for this fea- 
ture. Anisotropic pumping can produce strongly polarized 
masers even in the absence of a magnetic field (Western & 
Watson 1983); this mechanism is thought to explain the 
tangential polarization pattern seen the VLBA maps of 
SiO masers around evolved stars (Desmurs, et al. 2000). 

4.3. The v=0 Masers 

In contrast to the situation for the v=l mascrs, the SiO 
masers in the ground vibrational state almost certainly do 
fulfill the requirement that gil 3> R and T, hence their 
polarization angles should be perpendicular or parallel to 
the magnetic field. 

The v=0 maser emission region is extended with re- 
spect to our synthesised beam, with peak brightness tem- 
peratures of about 2000 K. To compute an upper limit 
for the stimulated emission rate, we assume that the 
masers radiate isotropically, so the radiation energy den- 
sity U = (4n/c)(2kT/X 2 ). Then R = BU = 1.4 x 10" 2 
s _1 . The spontaneous radiative rate A 2 i = 2.9 x 10~ 5 
s _1 is relatively slow, so the decay rate T from the 3=2 
level is set by the SiO-H2 collision rate. Modeling by Zcng, 
Lou, & Sun (1987) predicts that v=0 SiO masers originate 
in gas with density 10 7 cm~ 3 . If the gas kinetic temper- 
ature is ^300 K, the rate coefficient given by Bieniek & 
Green (1983) implies a collision rate of 6 x 10~ 3 s _1 . For a 
1 mG field the Zeeman frequency gfl/2n ~ 0.2 s _1 . Since 
g$l » R > T, the v=0 masers should be excellent tracers 
of the magnetic field direction. Reassuringly, the position 
angle of the SiO v=0 polarization appears to be stable with 
time. The maps in Figure were made from data taken in 
2000 December and 2001 February; the position angle is 
the same, P. A. ^80°, in lower quality images from 1996, 
1997, and 1999. We note that this is close to the median 
position angle of the v=l masers. 

One glaring discrepancy complicates the analysis: from 
observations made in 1994, Tsuboi et al. (1996) measured 
polarization position angles P.A. ~150° for the J=l-0 
masers in the v=0, 1, and 2 vibrational levels. As shown 
in Figure , these position angles differ by roughly 70° from 
the values we measure for the J=2-l line. Given that the 
position angle we measured for the v=0 J=2-l line was 
stable for 4 years, we doubt that this discrepancy can be 
attributed to time variability. It's also unlikely that the 
39" beam used for the J=l-0 measurements picked up 
a halo of emission polarized almost perpendicular to the 
core, since we find that virtually all of the polarized J=2~l 
flux originates in the central 3" region. This leaves two 
possibilities: either (1) the J=2-l and J=l-0 masers have 
intrinsically different polarizations; or (2) the intrinsic po- 
larizations are identical, but the observed position angles 
are twisted by Faraday rotation. Possibility (1) is hard to 
rule out. However, in the grid of theoretical calculations 
done by Nedoluha & Watson (1990) there appear to be no 
cases where the the J=2-l and J=l-0 polarization posi- 
tion angles differ by more than 20°. And, if the J=2 1 
position angle we measure is parallel or perpendicular to 
the magnetic field, then the field is aligned neither with 
the bipolar outflow from IRc2 (P.A. ^145°, as defined by 
the VLBA observations of the SiO masers), nor with the 
axis of the v=0 disk (P.A. -50°). 
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4.4. Faraday rotation 

We now explore the possibility that Faraday rotation 
causes the discrepancy between the J=2-l and J=l-0 po- 
larization position angles. Faraday rotation occurs when 
linearly polarized emission propagates through an ion- 
ized plasma with a magnetic field component along the 
line of sight. The plane of polarization is rotated by 
9 = X 2 RM, where the rotation measure is given by 
RM = 8.1 x 10 5 J n e B\\ dl rad m~ 2 . Here n e is the electron 
density in cm" 3 , Bn is the component of the magnetic field 
along the line of sight in gauss; and dl is in parsecs. If the 
propagation vector is along the magnetic field direction, 
then the plane of polarization rotates in a clockwise direc- 
tion, following the right hand rule. The minimum rotation 
measure which can bring the J=l-0 and J =2—1 position 
angles into agreement is RM = 3.3 x 10 4 rad m~ 2 , which 
rotates the J=2-l (A3. 5 mm) position angle by 23°, and 
the J=l (A7 mm) position angle by 92°. 

Because the Faraday rotation must be reasonably uni- 
form across the ~3" maser emission region, we presume 
that it occurs in the foreground Orion A H II region (M42), 
not in ionized gas local to IRc2 - e.g., on the surface of a 
circumstellar disk. What is the plausible rotation measure 
through Orion A? We estimate the electron column density 
from the work of Wilson & Jager (1987), who modeled the 
H II region as a series of 9 concentric cylindrical slabs in 
order to fit radio continuum and recombination line data. 
The line of sight to IRc2 passes through slabs 3-9, for 
which the model column density J n e dl ~ 650 cm -3 pc. 
We take the line of sight component of the magnetic field 
from the work of Troland, Heiles, & Goss (1989), who 
measured the 21 cm Zeeman effect in the sheet of neutral 
absorbing gas in front of Orion A. The inferred field is rea- 
sonably uniform, with no reversals of direction; the value 
toward IRc2 is approximately — 50/iG. The negative sign 
indicates that the field lines point toward us. Then the ro- 
tation measure is positive, meaning that the observed posi- 
tion angle increases with wavelength (Heiles 1987), which 
is what we are hoping to find. If the magnetic field is the 
same inside the H II region, then RM = 2.6 x 10 4 rad m~ 2 , 
within 20% of the value which brings the position angles 
into concordance. 

The estimated rotation measure is comparable with the 
the upper limit RM < 2.8 x 10 4 rad m~ 2 deduced by Rao, 
et al. (1998) from the agreement of the polarization po- 
sition angles for dust emission at wavelengths of 1.3 mm 
and 3.5 mm. We note that the majority of these dust po- 
larization detections were made west and north of IRc2, 
farther from the center of the H II region. 

4.5. Magnetic Field Direction 

Suppose we accept the hypothesis that Faraday rotation 
has twisted the J=2 1 position angle by +23°. As shown 
in Figure , this means that the intrinsic polarization angle 
is 57°, roughly parallel with the plane of the SiO "disk," 
and perpendicular to the axis of the high velocity outflow 
from source I. The inclination angle of the disk is uncer- 
tain, but the 3:1 axial ratio of both the SiO and H2O maser 
distributions suggest that its axis is tilted by < 20° from 
the plane of the sky, that is, that the disk is seen nearly 
edge-on. 

The magnetic field can be parallel or perpendicular to 



the SiO polarization vectors - but which is more likely? 
If the maser radiation travels more or less parallel to the 
magnetic field, then it should be polarized parallel to the 
projected field direction. Thus, our measurements are con- 
sistent with a magnetic field in the plane of the disk, as 
long as the field is oriented roughly along the line of sight. 
Obviously, such a field cannot be azimuthally symmetric 
(e.g., ring-shaped or spiral) around IRc2. If it were, the 
field would become perpendicular to the line of sight in the 
outer parts of the bowtie, and the polarization direction 
would flip by 90°. 

Our data are more easily reconciled with a field that is 
everywhere perpendicular to the plane of the disk. In this 
case the maser radiation propagates at an angle of 70° 
to the field, so the polarization should be perpendicular 
to the field. Then the field, at P. A. 145°, is parallel to 
the axis of the high velocity outflow and to the large scale 
Orion field mapped by far infrared observations (Schleun- 
ing 1998). The infrared measurements of Chrysostomou, 
et al. (1994) and Aitken et al. (1997) give a similar mag- 
netic field direction toward source I, along the SE edge of 
the IRc2 infrared peak. Our data provide no convincing 
evidence that the field has been pinched or twisted within 
~ 700 AU of source I. The twist seen in the infrared data 
occurs on slightly larger scales. 

5. CONCLUSIONS 

We attempted to infer the magnetic field morphology 
toward Orion-IRc2 from polarization measurements of the 
v=0 and v=l J =2-1 SiO masers associated with this star. 
The polarization vectors are expected to be perpendicular 
or parallel to the projected magnetic field direction pro- 
vided that the Zeeman splitting gfl 3> R and T, where R 
and T are the stimulated emission and decay rates from 
the maser's upper rotational level. 

We find that the polarization position angles of the v=l 
masers vary on time scales of years. These intense masers, 
which originate only ^40 AU from the star, probably are 
not suitable indicators of the magnetic field direction be- 
cause the stimulated emission rate exceeds the Zeeman 
splitting. Typically their fractional polarization is a few 
percent, although in 1999 a weak feature with 30% polar- 
ization appeared. 

The v=0 masers are much weaker and should be good 
probes of the magnetic field direction. These masers orig- 
inate 100-700 AU from the star in an elongated, bowtie 
shaped region, possibly along the top and bottom surfaces 
of an underlying flared disk. Maps of the v=0 emission 
with 0.5" (250 AU) angular resolution show that, contrary 
to our expectations: (1) the polarization vectors within 
~700 AU of the star are straight, with no evidence of a 
twist or hourglass morphology; (2) the measured polariza- 
tion position angle of 80° is neither parallel nor perpen- 
dicular to the high velocity outflow from IRc2; and (3) the 
J=2-l polarization position angle differs by 70° from the 
J=l-0 position angle previously measured by Tsuboi et al. 
(1996). 

The discrepancy between the J=2-l and J=l-0 position 
angles may be attributable to Faraday rotation by plasma 
in the foreground Orion A H II region. A rather high, but 
not implausible, rotation measure RM = 3.3 x 10 4 rad m~ 2 
is required to bring the two position angles into agreement. 
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Then the intrinsic polarization direction is parallel to the 
plane of the disk and perpendicular to the axis of the high 
velocity outflow from source I. Probably the polarization 
vectors are perpendicular to the field in this case, so the 
field threads the disk poloidally, parallel to the outflow 
and to the large scale magnetic field in Orion. The ab- 
sence of a discernable pinch or twist in the field near IRc2 
may indicate that as the molecular cloud collapsed the gas 
density in the pseudodisk exceeded 10 11 cm -3 , at which 
point the magnetic field decouples from the infalling gas 
(Nakano & Umebayashi 1986). Alternatively, the outflow 



or subsequent evolution of the disk has erased all traces of 
the pinched field. 

The high fractional polarizations observed for the v=0 
masers - up to 50% in our J=2-l data and up to 80% 
in the J=l-0 observations of Tsuboi et al. (1996) - sug- 
gest that anisotropic excitation by infrared photons play 
an essential role in pumping these masers. 

This work was partially supported by NSF grant AST- 
9981308 to the University of California. 
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Fig. 1. — Total intensity, fractional polarization, and position angle for the v=l J=2-l SiO masers (86.243 GHz) in Orion-IRc2, at 6 epochs. 
Fractional polarizations and position angles are plotted as points, with ±1(X error bars indicated. A dashed line is drawn at P. A. = 80°; this 
is the position angle of the v=0 masers. 
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Fig. 2. — Polarization vectors plotted at the fitted centroid positions for the v=l J=2-l SiO masers in Orion-IRc2, for 1 km s _1 velocity 
channels. Data for 3 epochs are overplotted to search for correlations between centroid positions and polarizations. Radio source I, denoted 
by a star, is located near the middle of the maser distribution. The arc of maser spots SE of source I corresponds to the blueshifted peak 
of the maser spectrum (Figure 1), while the arc of spots to the NW corresponds to the rcdshifted peak. Position uncertainties are of order 
±0.01". 
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Fig. 3. — Total intensity and polarized intensity v=0 J=2-l SiO spectra integrated over 30" X 30" and 3" X 3" boxes centered on Orion-IRc2. 
Spectra for the 30" box were generated from BIMA C array data taken in 2001 March, with a 10.4" X 6.7" synthesized beam; for the 3" 
box, from A and B array data taken in 2000 December and 2001 February, with a 0.97" X 0.50" beam. The polarized intensities (two lower 
curves) are similar, indicating that most of the polarized flux originates in the 3" inner region. 
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4. — Channel maps of J=2-l v=0 SiO emission in Orion-IRc2, 



RA offset (arcsec) 



RA offset (arcsec) 



Fig. 4. — Channel maps of J=2-l v=0 SiO emission in Orion-IRc2, generated from BIMA A and B array data obtained in 2000 December 
and 2001 February. The total intensity is indicated by contours, running from 0.5 Jy/beam (325 K) to 3.25 Jy/beam (2100 K) in steps of 
0.25 Jy/beam. Polarization vectors are overlaid wherever the polarized flux density P > 0.2 Jy/beam and, simultaneously, the total intensity 
/ > 0.5 Jy/beam. The maximum fractional polarization is ~50% in the —5 and +17.5 km s — 1 maps; it is much weaker at line center. Offsets 
are relative to radio source I at a = 5 h 35 m 14?505, <5 = — 5°22'30'/45 (J2000). The LSR velocity, in km s — 1 , is shown at the upper right of 
each panel, the 0.5" FWHM synthesized beam at the lower left. 
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Fig. 5. — Total intensities (thick curves), polarized intensities (thin curves) and polarization position angles for the J=2— 1 and J=l— SiO 
masers in the v=0 vibrational state, toward Orion-IRc2. The J=2-l transition was observed interferometrically with BIMA (this paper). To 
average over the maser emission region, we convolved channel maps of the I, Q, and U Stokes parameters with a 3" FWHM Gaussian, then 
computed the fluxes and position angles toward source I from these smoothed maps. The J=l— data (Tsuboi et al. 1996) were obtained with 
the Nobeyama 45-m telescope with a 39" beam. The position angles measured for the J=2— 1 and J=l-0 transitions differ by roughly 70°. 
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Fig. 6. — (a) Polarization map of v=0 J=2— 1 SiO emission in Orion-KL, integrated over a 28 km s — 1 wide velocity range centered at 
Vlsr 5 km s — 1 . The contour interval is 0.3 Jy/beam, or 200 K. Vectors indicate the polarization direction wherever there is a >4 a detection 
of the polarized flux, corresponding to a position angle uncertainty of ±7°. The highest fractional polarization is 28%. The position of radio 
source I is denoted by an open star at the center of the map. (b) Same, except the polarization vectors have been rotated by —23° to correct 
for Faraday rotation, as discussed in the text. Thin lines at P. A. —40° and 150° indicate the axes of the high velocity outflow from source I, 
as defined by the cones of v=l masers close to source I (Greenhill et al. 1998). 



